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The  detonation  velocity  of  solid  explosives  contained 

in  long" cylindrical  tubes  of  varying  weight, 

by  Dr.  H.  Jones. 


Introduction, 

The  Chapman-Jouguet  theory  of  the  detonation  of  gases  contained  in 
rigid  tubes  has  been  applied  by  many  authors  to  the  detonation  of  solid  explos¬ 
ives.  Calculations  by  A.  Schmidt,  R«  Becker,  Langweiler,  and  more  recently 
Kistiakowsky  and  Wilson,  have  used  the  theory  in  sin  inverted  sense  to  calculate 
fraa  the  observed  detonation  velocities  the  covolumes  of  the  gases  formed  in 
the  detonation  process.  In  the  earlier  calculations  of  A.  Schmidt  and  Lang- 
weiler  the  covolumes  so  determined  seated  to  bear  little  relation  to  the  actual 
size  of  the  molecules.  Kistiakowsky  and  Wilson,  however’,  using  a  more  flex¬ 
ible  form  for  the  equation  of  state,  obtain  values  which,  reduced  to  roan 
temperature,  are  not  incanpatible  with  the  b  of  van  der  Wash's  equation  as 
determined  at  ordinary  pressures.  A  practical  aim  of  these  calculations  is 
the  prediction  of  the  detonation  velocity  of  possible  new  explosives.  Applied 
to  explosives  which  detonate  with  the  development  of  very  high  temperatures 
these  calculations,  regarded  from  this  practical  standpoint,  have  considerable 
success;  they  fail  completely  when  applied  to  ammonium  nitrate  or  similar  ex¬ 
plosives,  always  predicting  (with  the  covolumes  determined  from  other  explos¬ 
ives)  a  velocity  far  in  excess  of  the  observed  value.  The  same  failure  occurs 
with  calculations  made  along  the  lines  of  R.C.  1 66  (cf.  §  5e  of  this  note.) 

It  seems,  at  first  sight,  remarkable  that  the  Chapman-Jouguet  theory 
in  no  way  involves  the  velocity  of  the  reactions  which  lead  to  the  breakdown 
of  the  explosive  molecule  and  the  fozmation  of  the  product  gases.  This  is 
because  the  theory  visualizes  a  steady  state  in  a  rigid  tube  and  connects  the 
properties  of  the  unexploded  gas  with  those  of  the  completely  reacted  gases. 

The  distance  which  separates  these  two  conditions  is  immaterial  to  the  theory. 
There  thus  appears  to  be  a  fundamental  distinction  between  the  process  of  de¬ 
tonation  in  gases  and  in  solids,  viz.  that  in  gases  the  confining  tube  remains 
intact  whereas  for  solid  high  explosives,  however  strong  the  confinement,  sane 
expansion,  if  not  complete  shattering,  always  occurs.  It  is  clear,  therefore, 
that  for  solids  the  distance,  in  the  steady  state,  which  separates  the  onset 
and  the  completion  of  the  reaction  must  play  a  part  in  determining  the  deton¬ 
ation  velocity.  It  has  been  shown  by  Boyes  that,  assuming  rigid  confinement, 
the  pressure  is  actually  higher  in  the  reaction  zone  than  at  the  point  where 
the  reaction  is  complete.  Hence,  if  the  length  of  the  reaction  zone  is  com¬ 
parable  with  the  diameter  of  the  charge,  an  appreciable  expansion  must  have 
occurred  before  the  whole  of  the  explosive  is  converted.  In  other  words,  part 
of  the  explosive  is  being  transformed  at  an  effectively  lower  loading  density 
and  therefore  the  velocity  may  be  expected  to  be  less  than  that  given  by  the 
Chapman-Jouguet  theory  assuming  infinitely  strong  confinement.  The  fact  that 
for  the  more  powerful  high  explosives  the  calculations  already  mentioned  are 
not  unsuccessful  must  be  due  to  the  very  high  rate  of  the  reactions  taking 
place,  so  that  in  general  the  length  of  the  reaction  zone  is  very  small  compared 
with  the  diameter  of  the  charge.  Nevertheless  reasons  are  given  in  §  5  for 
believing  that  in  T.  N.  T.  the  observed  velocities  may  be  as  much  as  10  per  cent, 
less  than  the  theoretical  values  assuming  rigid  confinement. 

In  a  recent  communication  (r.C.193)  Professor  G.  I.  Taylor  has 
analysed  the  process  of  an  explosion  in  a  long  cylindrical  tube  end  has  shown 
how  the  radius  of  the  tube  increeises  under  the  adiabatic  expension  of  the  pro¬ 
duct  gases  of  the  detonation  wave.  The  analysis  of  the  present  note  is  along 
the  lines  of  Professor  Taylor's  paper;  the  main  difference,  in  the  mathematics, 
is  that  in  Professor  Taylor's  calculation  there  existed  a  known  relation  between 
the  pressure  and  the  density,  viz.  the  adiabatic  relation,  whereas  in  the 
present  calculations  such  an  algebraic  relation  does  not  initially  exist  but  is 
replaced  by  a  differential  relation.  In  other  words,  in  the  present  calcul¬ 
ations,  the  differential  equation  expressing  the  adiabatic  condition  cannot  be 
integrated  independently  but  only  simultaneously  with  the  equation  of  motion. 


g  1.  The  general  equations  of  the  detonation  process. 

We  consider  the  detonation  of  a  long  cylindrical  charge  surrounded 
by  a  casing  of  weight  X  per  unit  length.  When  the  detonation  process  is 
fully  developed  a  steady  state  is  established.  If  the  detonation  velocity 
is  U  we  consider  the  process  fran  axes  moving  with  this  velocity.  The  shock 
wave  front  in  the  solid  is  then  at  rest.  This  point  is  taken  as  the  origin 
x  s  0.  Behind  the  shock  wave  front  the  velocity  of  the  reacting  material 
we  denote  by  q  which  is,  of  course,  a  function  of  x.  There  is  a  discontinu¬ 
ous  change  of  pressure  and  density  at  the  shock  wave  front.  Let  Eu(p,f>)  be 
the  energy  per  unit  mass  of  the  unexploded  material  at  pressure  p  and  density 
reckoned  frcm  a  zero  at  p  =  0.  The  equations  connecting  the  conditions 
on  each  side  of  the  discontinuity  are  then 

r  r-q°  *  UA 

•  -  UA  (U  -  q.)  =  uaA  (1  -A)  .  (1) 

.  *.<w.)  •  Ki-jt)  - 

where  A  is  the  unexploded  density,  and  f0,p0  ,  q<,  are  the  density,  pressure, 
and  velocity  in  the  shock  wave  front. 

The  equation  of  continuity,  when  the  angle  of  the  expanding  tube 
is  small,  can  be  expressed 


Raq f>  m  R^JA  .  (2) 

when  R„  is  the  initial  tube  radius  and  R  the  radius  at  the  point  x.  The 
equation  of  motion  to  the  same  approximation  is 

q  sf  *  0  .  (3> 

Let  Ee(p,p)  denote  the  internal  energy  per  unit  mass  of  the  converted 
explosive  material,  and  let  n  denote  the  fraction  of  the  explosive  which  has 
been  converted  at  a  particular  point  x.  Then  if  E  be  the  total  energy  per 
unit  mass 


E  c  nEe  +  (1  -  n)Ea 


(4) 


The  pressure  in  both  parts  of  E  mast  be  the  same  but  the  density  may,  of 
course,  be  different. 

Also  we  denote  by  H0  the  whole  of  the  chemical  energy  released  per 
unit  mass  when  the  reaction  is  complete.  At  any  point  x  the  energy  released 
is  denoted  by  H  and  therefore 

H  =  nH0  .  (5) 


In  addition  to  the  kinematics!  equation  of  continuity  (2)  and  the 
dynamical  equation  (3)  we  hare  also  a  thermodynamical  equation  expressing  the 
conservation  of  energy.  We  write  this  in  the  form 

X(E_H).p£(i)  .  0  .  <6> 

This  equation  assumes  that  no  heat  is  transferred  frcm  one  element  to  the  next 
during  the  development  of  the  reaction  frcm  the  shock  ware  front  to  the  point 
at  which  the  reaction  is  complete.  We  are  therefore  emitting  any  consider¬ 
ation  of  the  actual  kinetics  of  the  process  by  which  the  reaction  is  propagated 
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with  the  result  that  the  problem  of  the  propagation  of  detonation  is  insol¬ 
uble  without  a  further  basic  assumption,,  In  the  case  of  gases  the  assumption 
generally  known  as  the  Chapnan-J ougue t  condition  has  proved  successful  and  has 
been  Justified  theoretically  by  Bqyesr  for  the  case  in  which  the  mechanism  of 
the  reaction  depends  on  heat  conduction. 

In  the  present  calculations,  which  attempt  to  take  account  of  the 
finite  time  required  for  the  completion  of  the  reaction  in  the  detonation 
wave  and  the  consequent  expansion  during  this  reaction,  we  apply  the  Chapnan- 
Jouguet  condition  in  a  slightly  generalised  form.  The  assumption  made  is 
that  at  the  point  where  the  reaction  is  complete  the  velocity  is  equal  to  the 
local  velocity  of  sound.  If  we  denote  by  X  the  value  of  x  at  the  point 
where  the  reaction  is  complete,  and  the  conditions  at  this  point  by  the  sub¬ 
script  i  the  above  assumption  is  expressed  by 

<1,  =  C,  . .  (7) 

where  C,  is  the  velocity  of  sound  at  x  -  X.  For  infinitely  strong  confine¬ 
ment  this  is  identical  with  the  usual  condition  applied  to  gases  and  solids. 
Since 

*  (pi  . « 

the  condition  (7)  becomes 


In  addition  to  equation  (3)  there  is  another  dynamical,  equation  determining 
the  expansion  of  the  casing,  viz. 


tR 

77*  " 


(10) 


in  which  t  has  been  replaced  by  x/U  since  we  are  dealing  with  a  steady  state. 


§  Reduction  of  the  equations  to  a  form  suitable  for  integration. 

In  order  to  solve  the  above  equations  we  make  a  number  of  assump¬ 
tions.  These  are  of  a  different  kind  from  (7);  they  express  the  initial 
chemical  information  upon  which  any  theory  of  detonation  must  be  based,  and 
could  be  modified  without  affecting  the  main  conclusions. 

Fir3t  we  assume  that  the  reaction  velocity  is  constant  throughout 
the  reaction  zone,  i.e. 


n  =  (x/X) 


(11) 


Secondly  we  take  a  very  simple  form  for  Ee,  viz. 

E*  =  c  £  .  (12) 

where  c  is  a  constant  of  the  nature  of  a  specific  heat.  Thirdly  we  assume 
that  H0  is  independent  of  p  and  p  (this  has  been  shown  to  be  very  nearly  the 
case  for  T.  N.  T.  ) 


x 

X 


H, 


H 


(13) 
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Fourthly  we  assume  that  in  the  expression  for  the  energy 
E  =  £c£+/l-*)E 

x  r  \  x/  u 

Eu  can  be  regarded  as  a  constant. 

It  is  convenient  to  use  the  following  abbreviations : - 

y  =  &  ,  h  =  2°,  and  r  =  ^/R0* 

C  0  uA 

Equation  (6)  can  then  be  written,  making  use  of  (l). 

Equation  (3)  with  (2)  becomes 


(14) 


(15) 


y 


*  £ 


(r  y  ) 


=  0 


which  can  be  rewritten 


(16) 


(17) 


and  integrating  from  0  to  x  and  using  equations  (1 )  we  have 


^  =  i  -  y  ♦  Mi  -  *•'")  +  ?  </> 


(18) 


where 


<f> 


X 


(19) 


It  is  convenient  to  write  the  equations  in  this  way  since  1  -  r 
is  small  for  heavy  confinement,  vanishing  of  course  for  rigid  tubes.  If  we 
integrate  in  a  similar  way  equation  (15)  we  obtain 


cxy  [  1  -  y  +  Ml"**'4)  +  ?  +  xjy  -  iya-  y0+  Mj  .  (2°) 

+  i  J  ( 1  -r^)y  dx  +  ^  dx  =  £  h0+  %( 1  -y0  )  J  * 

A  partial  integration  gives 

j^||dx  =  (f>y  -J  y(l-r*)  ||  dx 
and  hence  (20)  becomes 

y3^  (X  +  2cx)  -  cx(l-r^)J  -  (2  +^)(cx  +  X)  y  .  (21 ) 

+  [  (X  +  2ht)x)  -  (X  -  x) ( 1  -  y 0)*J  = 


0 
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This  equation  is  not  merely  a  quadratic  in  y  since  i>  also  depends  on  y.  It 
is,  however,  a  useful  form  for  our  purpose  since  fi  is  small  for  heavy  confine¬ 


ment,  i.e„  large  X  • 

Fran  equation  (8)  for  the  velocity  of  sound  we  find 

(«/*  ufeOjr  .  («> 

and  thus  the  Chapman -J ouguet  condition  (9)  with  equation  (2)  gives 

(s'/-  (*/•  A"  -«%(¥> . («) 

which  earn  be  rewritten  with  the  help  of  (18) 

y,  [(l+2c)  -  0(1-0  j  =  (l+c)[  1  +  fa,  (l-r,"4)+  .  (24) 

when  all  quantities  refer  to  the  point  where  the  reaction  is  complete.  Also 
at  x  =  X  (21)  becomes 

y*  [  (l+2c)  -  c(l -r"*)J  -  (2+^)(l+c)y,  +  (l+2hj  -  0  .  (25) 

which  with  (24)  gives 


2he 


25 


ca  »  (1+c)*  [  (f,  +  j  p*-  jy/(l-Oa]  ♦  c(l-r,~4) 
(l+2c)  -  o(1-r,**) 


(26) 


For  infinitely  strong  confinement 


^  =  (l-r_lf)  =  0  and  we  have 


2H0 


(27) 


where  U*,  denotes  the  detonation  velocity  in  this  limiting  case. 
We  can  write  equation  (26)  therefore  as  follows 


This  equation  does  not,  of  course,  give  U  until  y  and  r  have  been  obtained  as 
functions  of  x.  This  is  achieved  by  successive  approximation  for  the  case  of 
heavy  confinement. 


§  3*  The  solution  for  a  non-exp andlng  tube. 

In  this  special  case  if)  =  1  -r^  =  0  and  2h0  =  ca/0+2c)  and  therefore 
(21 )  reduces  to 


y*(1  ♦  2cf)-2(cf.1)y.[l  -  <1-f  )(1-y.)‘J  =0 


(29) 
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We  denote  the  solution  of  this  equation  by  y  °  (x)  and  regard  it  as 
the  zero  order  approximation  for  the  solution  in  the  case  of  an  expanding  tube. 
The  solution  of  (29)  which  satisfies  the  initial  condition  that  y  =  y  for 
x  =  0  is 


y(o)(x)  = 


1  +  c 


4 


2cf  + 


(30) 


When  x  =  X  we  have  y,°^  a  (l+o)/(2c+1 )  which  is  independent  of  yD.  In  order  to 
see  the  nature  of  the  variation  of  yw  with  x  we  may  take  a  special  example 
c  =  ?  which  corresponds  to  y,(o)  =  %  or  ( P ,®/  A  )  =  1  *33*  Instead  of  y  we 
give  1-y  which  is  equal  to  p  /uaA. 


V\(a/X) 

(a/^)\^ 

0 

.2 

•  4 

.6 

.8 

1.0 

.75 

.25 

.25 

.25 

.25 

.25 

.25 

060 

.40 

.388 

.374 

.356 

.330 

.25 

•  50 

.50 

*  474 

.  446 

.413 

.369 

.25 

.40 

.60 

•  563 

.516 

.467 

.406 

.25 

It  will  be  seen  that  the  pressure  is  always  greater  in  the  reaction 
zone  than  at  the  point  where  the  reaction  is  complete.  There  is  no  solution 
possible  when  (&  / f>0  )  >  (l+c)/(l+2c),  that  is  it  is  necessary  that  po  f>{ 
in  order  that  a  solution  may  exist.  When  f>0  =  the  density  and  the  pressure 
are  constant  throughout  the  reaction  zone. 


(o| 

The  rate  of  change  of  y  with  x  at  the  point  x  = 
following  equation  obtained  from  (30) 


X 


s 


( )  (i£ [  - 

2(i-y0) 


0  is  given  by  the 


(31) 


In  the  following  we 


denote  X 


by  a. 


§  4.  An  approximate  solution  for  an  expanding  tube  when  the  expansion  is  small. 

We  obtain  an  approximate  solution  to  (21)  by  calculating  (l-r"4)  and 
r  from  the  known  zero  approximation  for  y,  viz.  y^(x).  The  first  approxi¬ 
mation  therefore  to  r  as  a  function  of  x  is  given  by  (10)  and  (18)  neglecting 
in  these  equations  the  terns  (l-r'^)  and  <p  .  Hence  the  equation  giving  the 
first  approximation  to  r  is 


2TTA 


(i  -  y  Co)(x)) 


(32) 


We  make  the  further  approximation,  which  the  preceding  table  shows  to  be  not  too 
unsatisfactory,  of  replacing  y^(x)  by  y0  +  a—.  The  solution  of 

X 


■  (l 


yc 


r 


^r 


•  ••••• 


(33) 
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to  first  powers  of  A/X  is 


r  =  1  ♦  J2A.  r 

\ 


!(fl 


and  therefore 


1  -  r  *  = 


and 


* 


4TT  A  X' 

A 


O-^dM  (il 


4T  A  *!  a 

3A 


Mur-HD 


(34) 


(35) 


(36) 


With  equs.  (35)  and  (36)  a  first  approximation  to  (21)  could  be  ob¬ 
tained.  Clearly  higher  approximations  could  only  be  found  numerically,  i.  e. 
for  a  particular  case  in  which  the  numerical  values  of  all  the  coefficients 
are  known.  It  would  then  be  possible  to  find  solutions  of  any  degree  of 

accuracy.  For  the  present  purpose  it  is  more  interesting  to  obtain  general 

results  for  the  case  of  heavy  confinement,  i.e.  for  small  1-r  tt.  To  first 
powers  of  h  Xa  therefore  equation  (28)  becomes 


1  ♦ 


(jf  2/{* 


(37) 


and  using  equations  (35)  and-  (36)  we  find 

=  1  +  .  (38) 


where 

f  M 


and  a  is  given  by  equation  (31).  Since  (/°/Va)  =  (2c,+1  )/(c+1 ),  the  factor 
f  can  be  expressed  in  terns  of  (/£  /A  )  and  (^°/A  ).  The  following  table  gives 
an  idea  of  the  magnitude  of  this  factor. 


If  t  denotes  the  time  taken  for  the  whole  conversion  of  the  explos¬ 
ive  then  X  =  U  f ,  and  if  also  in  place  of  A  we  use  the  ratio  of  the  weight 
of  the  case  to  the  weight  of  the  explosive  <*-  so  that  X  =  TT  R*  , 
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equation  (38)  may  be  written 


(40) 


Still  another  way  in  which  the  result  can  be  shown  is  as  follows 


u<o  -u  _  f 

u«  4  R„  <r 


(41) 


where  cr  is  the  mass  per  unit  area  of  the  encasing  tube.  This  would  be  applic¬ 
able  to  the  case  of  the  variation  of  U  with  diameter  for  tubes  of  given  thick¬ 
ness.  The  relation  between  «•  and  R0  shewn  in  (40)  has  already  been  noted  by 
Professor  Taylor  (R.  C.  193). 


§  5-  Comparison  with  experimental  data. 

There  does  not  appear  to  be,  unfortunately,  much  published  data  in 
which  sufficient  details  are  given  to  make  a  quantitative  comparison  with  the 
foregoing  calculations.  There  are,  however,  a  number  of  properties  of 
detonating  explosives  which  appear  to  fit  in  fairly  well  with  the  above 
equations. 

(a)  It  is  observed  that  although  the  detonation  velocities  vary  with 
diameter  and  degree  of  confinement  the  variations  are,  in  the  case  of  the  more 
powerful  explosives,  only  of  the  order  of  a  few  per  cent,  for  comparatively 
wide  variations  of  oc  and  Rq.  It  follows,  therefore,  according  to  (40),  that 
X  =  Ut  must  be  appreciably  less  than  the  diameter  of  the  charges  ordinarily 
used,  e.g.  it  must  not  be  more  than  1  or  2  mm.  for  substances  such  as  picric 
acid  or  T.  N.  T. 

It  is  very  tempting  to  correlate  the  excess  width  of  the  well  marked 
trace  on  moving  film  photographs  with  the  thickness  of  the  reaction  zone.  A 
reason  supporting  this  hypothesis  is  that  the  comparative  sharpness  of  the 
trace  makes  it  very  difficult  to  assign  it  to  temperature  radiation  since  the 
gases  for  same  distance  behind  the  detonation  wave  front  have  still  very  high 
temperatures.  If  not  to  temperature  radiation  it  would  appear  to  be  due  to 
chemi-luminescence,  i.  e.  to  the  radiation  emitted  during  the  breakdown  of  the 
explosive  molecule  and  the  formation  of  the  gas  molecules.  P.  Laffitte 
(Ann.  de  Physique,  1925,  10,  p«659)  remarks  that  in  all  his  photographs  the 
trace  was  quite  sharp  and  that  the  thickness  of  the  image  corresponded  closely 
with  the  thickness  of  the  tube  photographed.  That  is  to  say  it  was  not 
appreciably  thicker  than  would  be  given  ty  an  infinitely  thin  disc  of  light 
of  the  diameter  of  the  charge.  If  we  accept  the  hypothesis  that  the  thickness 
of  the  disc  of  light  which  is  photographed  in  these  experiments  can  be  identi¬ 
fied  with  the  width  of  the  reaction  zone  X,  Laffitte 's  experiments  show  that 
X  <<  Re,  which  is  necessary  according  to  (40)  if  the  detonation  velocity  is  not 
to  be  far  too  dependent  on  R0  and  cx.  . 

(b)  Such  quantitative  experimental  work  a3  has  been  published  on  the 
variation  of  detonation  velocity  with  diameter  of  charge  refers  mainly  to 
charges  with  very  light  paper  confinement.  The  effect  of  the  inertia  of  the 
air  shock  wave  is  much  more  difficult  to  consider  than  that  of  a  metal  casing. 
However,  if  we  replace  cr  by  an  effective  surface  density  cr«^  due  to  the 
outside  gas  and  the  light  covering,  it  makes  a  very  rough  comparison  with  (4-1 ) 
possible. 


R.  FOrg  (S.S. ,  1916,  p.37)  has  given  the  following  figures  for  the 
variation  of  detonation  velocity  for  T.  N. T.  enclosed  in  paper  covers. 
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Diameter 

Cubic  density 

Det.  velocity 

(iW 

2.0  cms. 

0.80  gn/cm3 

3765  ny'sec. 

.OBI 

3.0  " 

0.80  " 

3905  " 

.047 

4.0  " 

0.79  M 

4054  " 

oOfl 

5.0  " 

0.83  * 

4099  " 

.00 

In  the  last  column  denotes  the  detonation  velocity  at  the  greatest  dia¬ 

meter  used,  viz*  5  «30S 

If  in  (41)  we  take  the  value  of  f/4  *  U  aT A  A  /cr^.  to  be  0.1  cms. , 
the  values  of  (U®  -  U)/U»  would  be  for  the  diameters  of  tne  above  table  0.1, 
.067,  .05,  .04  which  gives  the  same  sort  of  variation  since  <  U®,  .  We 
have  seen,  however,  that  f  is  of  the  order  of  4  and  for  the  above  case 
A  =  .8  gp/rm3  and  if  UT:  is  taken  to  be  not  more  than  1  ran. ,  as  photographs 
suggest,  then  cr^  is  not  more  than  0.1  gn/cm*,  which  is  a  plausible  value 
being  not  much  more  than  the  weight  of  a  paper  covering. 

The  effect  of  the  variation  of  the  diameter  of  charge  on  the  deton¬ 
ation  velocity  is  still  more  pronounced  for  the  less  powerful  explosives  such 
as  ammonium  nitrate  or  potassium  nitrate  mixtures  or  sane  of  the  'milder' 
dynamites.  For  example,  Dautriche  found  for  cheddite  an  increase  of  about 
2C ft£  in  the  detonation  velocity  when  he  increased  the  diameter  fran  2  to  4  cans. 
Also  Laffitte  finds  a  detonation  velocity  for  No. 1  dynamite  about  500  m/sec. 
less  than  the  value  found  by  Dautriche  for  the  same  substance  and  ascribes  the 
difference  to  the  fact  that  his  tubes  were  4  mm.  diameter  whilst  Dautriche 
used  20  nm.  tubes.  The  increased  effect  of  the  diameter  or  the  casing  weight 
for  the  weaker  explosives  is  in  agreement  with  equation  (41 )»  since  for  these 
explosives  the  temperature  of  the  detonation  wave  is  lower  than  for  the  power¬ 
ful  high  explosives  and  consequently  T  is  larger. 


(c)  Calculations  of  the  detonation  velocity  of  T.N. T.  based  on  Bridgpan‘s 
compressibility  data  for  nitrogen  (R.C. 1 66)  were  found  to  be  from  12  to  15  per 
cent,  greater  than  the  observed  values.  These  calculations  assumed  infinitely 
strong  confinement  and  it  appears  that  part  at  least  of  the  discrepancy  may 
be  due  to  this  assumption,  At  a  loading  density  of  0.8  gn/cm'?  the  above  dis¬ 
cussion  suggests  that  U  might  be  10  per  cent,  less  than  Uoo  for  charges  not 
exceeding  1  cm.  diameter  and  with  very  light  covering. 


(d)  A  characteristic  feature  of  detonation  velocity  -  loading  density 
curves  is  a  tendency  to  flatten  at  high  loading  densities,  i.e.  for  the  rate 
of  increase  of  the  detonation  velocity  with  density  to  became  smaller.  In 
the  case  of  some  explosives,  notably  cheddite,  the  curve  actually  reaches  a 
maximum  beyond  which  the  velocity  decreases  as  density  increases.  It  was 
shown  in  R.C.  166  that,  in  general,  the  temperature  of  the  detonation  wave 
decreased  with  increasing  loading  density.  Although  the  reaction  mechanism 
of  the  detonation  wave  is  not  understood  at  present  it  is  a  fair  assumption 
that  the  rate  of  reaction  will  decrease  with  decreasing  temperature;  in  other 
words  "t  will  increase  as  the  temperature  falls.  This  could  therefore  explain 
the  flattening  of  the  U,A  curve  which  does  not  appear  to  be  due  to  the 
equation  of  state,  for  in  none  of  the  cases  so  far  examined  has  the  U®  , A 
curve  shown  any  such  tendency  at  high  densities. 

A  limitation  to  the  applicability  of  (38)  (apart  from  the  many 
approximations  already  discussed)  is  that  X  ,  in  the  way  it  has  been  used  in 
these  calculations,  cannot  in  practice  be  increased  indefinitely.  Merely 
increasing  the  weight  and  thickness  of  the  confining  tube  does  not,  beyond  a 
certain  point,  increase  A  because  of  the  compressibility  and  fluidity  of  all 
actual  substances.  Thus  the  fact  that  the  detonation  velocity  of  T.N.T.  at 
high  densities  does  not  appear,  judging  by  published  data,  to  exceed  7,000 
m/aec.  even  with  fairly  strong  confinement  does  not  necessarily  mean  that  this 
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velocity  is  the  ideal  maximum  corresponding  to  infinitely  strong  tubes,  and 
does  not  therefore  rule  out  the  suggested  explanation  of  the  observed  dimin¬ 
ishing  of  (dU/d  A  )  for  high  A  . 

(e)  Using  the  methods  of  R.C.  166  the  ideal  maximum  detonation  velocity 
Uoo  of  ammonium  nitrate  as  a  function  of  loading  density  has  been  calculated. 
The  results  axe  given  below. 


A  gycxa3 

0.475 

0.  580 

0.708 

0. 837 

0.984 

U^  m/sec. 

2608 

3069 

3781 

4570 

5450 

T°K» 

1431 

1430 

1425 

1400 

1300 

It  is  interesting  that  observations  of  detonation  velocities  of 
ammonium  nitrate,  and  of  ammonium  nitrate  with  small  quantities  of  other  sub¬ 
stances,  vary  very  extensively  between  different  observers  and  appear  to  be 
very  sensitive  to  the  conditions  of  confinement,  and  even  more  to  the  addition 
of  small  quantities  of  oxidisable  material,  e.  g.  carbon  or  aluminium.  The 
actual  values  of  U  (cf.  Marshall,  Explosives)  are  considerably  less  than  the 
values  of  U®  given  above  at  the  corresponding  loading  density.  These  facts 
appear  to  fit  nicely  with  the  very  low  temperature  of  the  detonation  as  shown 
in  the  above  table  which  indicates  an  exceptionally  large  reaction  time  and 
therefore,  according  to  (41)  an  exceptional  sensitivity  to  the  conditions  of 
confinement.  The  addition  of  carbon  or  aluminium  allows  for  the  absorption 
of  the  excess  oxygen  in  the  detonation  products  of  ammonium  nitrate  and  the 
formation  of  CO*  or  AI2O3.  The  great  heats  of  formation  of  these  molecules 
result  in  a  considerable  increase  in  the  temperature  of  the  detonation  wave 
and  therefore  to  a  reduction  in  the  reaction  time,  which  according  to  (41 ) 
brings  the  detonation  velocity  closer  to  U«  . 
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